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We study the coherent flow of a guided Bose-Einstein condensate incident over a disordered 
region of length L. We introduce a model of disordered potential that originates from magnetic 
fluctuations inherent to microfabricated guides. This model allows for analytical and numerical 
studies of realistic transport experiments. The repulsive interaction among the condensate atoms 
in the beam induces different transport regimes. Below some critical interaction (or for sufficiently 
small L) a stationary flow is observed. In this regime, the transmission decreases exponentially with 
L. For strong interaction (or large L), the system displays a transition towards a time dependent 
flow with an algebraic decay of the time averaged transmission. 
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I. INTRODUCTION 

The extraordinary experimental control achieved over 
atomic Bose-Einstein condensates (BEC) provides new 
testgrounds for phenomena issued from many different 
fields. On the one hand these systems allow to study ex- 
tensively nonlinear phenomena such as four- wave mixing 
, propagation of bright @] and dark 0| solitons or the 
dynamics of Bloch oscillations in presence of atom-atom 
interactions 0, [B| . On the other hand the rapid progress 
in this field has lead to a number of fascinating exper- 
iments probing complex condensed matter phenomena, 
such as the Mott transition in optical lattices [(J, the 
creation of vortices 7], the Josephson effect Q or the 
BEC-BCS crossover [jj. Bose-Einstein condensates link 
these two prominent fields of current research in a excit- 
ing and unique way. 

Wave mechanical transport in atomic vapors appears 
as a new direction for these trans-disciplinary studies 
that provide deeper insights into transport phenomena in 
presence of interaction. Indeed, BEC systems are intrin- 
sically phase coherent, as are the clean two-dimensional 
electronic structures studied in mesoscopic physics at low 
temperatures. Besides, interaction is much more simply 
modeled in BEC systems than the electrostatic electron- 
electron potential and its sign (repulsive or attractive) 
and strength can be tuned almost at will. The link be- 
tween matter-wave physics and electronic transport phe- 
nomena became ultimately apparent with the advent of 
microscopic traps and waveguides for the atoms, known 
as atom chips [lfj, [llj, [12( • Related studies include the 
attempt to generalize Landauer's theory of conductance 
to cold atoms |l3j . the atom blockade phenomenon in 
quantum-dot like potentials 14] as well as nonlinear reso- 
nant transport of Bose Einstein condensates , to men- 
tion just a few examples. 

A new direction in this context is the transport of Bose- 
Einstein condensates through disordered potentials. A 
relevant question is to which extent a Bose-Einstein con- 



densate is subject to Anderson localization jla. Il7| in 
presence of disorder, as well as how this scenario is af- 
fected by the atom-atom interaction. There is a grow- 
ing interest in the BEC community for issues related to 
the behavior of matter waves in disordered potentials. 
It started with the observation of "fragmentation of the 
condensate" over a microchip [l^. Nowadays a random 
potential is routinely engineered using an optical speckle 
pattern and its effects over the expansion of the conden- 
sate have been explored in Refs. pH |2(1I]. 

In contrast to studies where the condensate is initially 
at rest, we focus in the present paper on the effect of 
disorder on a propagating Bose-Einstein condensate. In 
an adiabatic approximation, the dynamics reduces to an 
effective one dimensional (ID) transport problem, this is 
the so called "ID mean field regime" |2l|. We further- 
more assume that the mean kinetic energy of the atoms in 
the condensate is larger than the typical height of the bar- 
riers induced by the disorder potential, i.e., perfect trans- 
mission is expected by classical mechanics. For the sake 
of concreteness, we restrict ourselves to one specific type 
of disorder: the one experienced by a condensate that is 
magnetically trapped above a corrugated microchip. To 
this end we introduce a model that could be character- 
ized as a "dirty wire model" where the current in the 
microfabricated wire has white noise fluctuations. This 
simple model captures most of the caracteristics of the 
random potentials observed over corrugated microchips. 
We point out, however, that our results are not expected 
to be sensitive to the particular type of disorder, as long 
as the latter is sufficiently smooth and can be character- 
ized by a well-defined correlation length. 

Previous theoretical studies of the effect of disorder on 
the transmission of nonlinear waves mainly focussed on 
attractive interaction and looked for stationary solutions 
of the problem (see the review [2^| ). Then, one has to 
choose between fixed input and fixed output boundary 
conditions. The latter case is less realistic, but simpler 
to discuss: It leads to algebraic decay of the transmission 
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|23j. The former case is complicated by the advent of 
multistability. However, the results of Knapp et al. [24[ 
show that, for short sample size, the mean transmission 
is poorly affected by a weak nonlinearity (as compared to 
the linear case), whereas for larger samples and stronger 
nonlinearity, evidences of derealization are found. . 

Realistic transport processes of Bose-Einstein conden- 
sates are different from the above mentioned studies be- 
cause they typically involve particles experiencing repul- 
sive interactions. We will see below that in this case the 
assumption of stationarity is not appropriate because for 
large disordered region or strong nonlinearity stationary 
solutions are dynamically unstable. In typical experi- 
ments the population of a given final state can only be 
achieved through a time-dependent process (such as the 
gradual filling of an initially empty waveguide with mat- 
ter waves). As a result, if a stationary scattering state 
is unstable, the transport properties of the condensate 
may be unrelated to the transmission coefficient associ- 
ated with that state, whereas a study of stationary flows 
might misleadingly give some weigth to this state (if the 
transmission is averaged over all possible stationary so- 
lutions for instance). 

We thus consider a setup that is relevant to experi- 
mental realizations and adapted to this specific trans- 
port scenario explicitely taking into account the possi- 
bility of time-dependent scattering: a coherent source of 
atoms emits matter waves that propagate in the mag- 
netic waveguide and encounter on their path a disorder 
region of length L. We show that the presence of a repul- 
sive atom-atom interaction has dramatic effects on the 
transport properties of the condensate. As is the case 
for attractive nonlinearity, Anderson localization is ob- 
served only in the regime of small interaction strengths 
and sample lengths. In this regime the transmission 
deacreases exponentially with L (oc exp(— L/Lj oc )), with 
a localization length Li oc modified by the interaction. 
For large sample lengths or strong interaction, time- 
dependent scattering processes occur. In contrast to the 
previous regime, one observes an Ohmic decrease of the 
time-averaged transmission (oc L^ 1 ). 

The paper is organized as follows. In Sec. [H] we set 
up the theoretical framework that is necessary to study 
transport through mesoscopic waveguides, introduce an 
effective one-dimensional Gross-Pitaevskii equation and 
present a numerical method that is particularly suited to 
study transport processes of Bose-Einstein condensates 
in waveguides. In Sec. lllll we introduce a one-dimensional 
model for the random magnetic potential along the center 
of the waveguide. We will show that a microscopic mean- 
dering of the current in the wire on the atom chip leads 
to a Lorentzian correlated random potential. In Sec. IIVI 
we investigate the regime of weak disorder potentials and 
give a simple analytic expression for the condensate wave- 
function in the guide. In Sec. |V] we discuss numerical 
results for transport through moderate and strong disor- 
der regions. We consider in particular the scaling of the 
transmission with the length of the disorder region. The 



paper closes with some concluding remarks. Some tech- 
nical points are given in the appendices. In Appendix 1X1 
we derive a relation between the mean transmission and 
the correlation function of the disorder potential. In Ap- 
pendix [B] we rederive, using standard WKB techniques, 
a result which is obtained heuristically in the main text. 

II. TRANSMISSION THROUGH WAVEGUIDES 

We consider a coherent beam of Bose-Einstein con- 
densed atoms at zero temperature, propagating through 
a cylindrical magnetic waveguide of axis x. The conden- 
sate is formed by atoms of mass m which interact via a 
two-body potential characterized by its 3D s-wave scat- 
tering length a sc . We consider the case of repulsive effec- 
tive interaction, i.e., a sc > 0. The condensate is confined 
in the transverse direction by an harmonic potential of 
pulsation ujj_. This transverse confinement is character- 
ized by the harmonic oscillator length aj_ = (h/muj^) 1 ^ 2 . 

In the following we restrict ourself to the "ID mean 
field regime" |2l| corresponding to a density range such 
that (a sc /a±) 2 <C nma sc <C 1, where niD denotes a 
typical order of magnitude of the ID density n(x, t) of 
the system. The first of these inequalities ensures that 
the system does not get in the Tonks- Girardeau limit 
and the second that the transverse wave function is the 
ground state of the linear transverse Hamiltonian, see, 
e.g., the discussion in Refs. [2ll |25|. In this regime the 
system is described by a ID order parameter ipix, t) (such 
that n(x, t) = \tp(x, t)\ 2 ) depending only on the spatial 
variable x along the guide. il){x : t) obeys the ID Gross- 
Pitaevskii equation 

dib ( h 2 d 2 \ 
lh ^={-^d^ 2+V{x)+anM )^ (1) 

with g = 2hw±a sc [2(||23,|23. V(x) is an effective one- 
dimensional potential along the waveguide, to which the 
condensate is exposed during the propagation process. 
We will see in Sec. EU how it may originate from irregu- 
larities of a wire used for creating the magnetic confine- 
ment. 

In absence of a potential (V(x) = 0) the plane wave 

ip(x, t) = ^riQ exp(ikx — i^it/K) (2) 

is obviously a solution of the Gross-Pitaevskii equation 
• It satisfies the dispersion relation 



where the particle current is given by J ~ nohk/m. 
Therefore, the chemical potential /i and the equilibrium 
constant density no of a freely propagating condensate 
beam are determined by the current </, the wave vec- 
tor k, and the effective interaction strength g. At this 
point, we mention that it was demonstrated in |28j | that 



3 



~- — 
> 

o 
-l 



Upstream region 


A Downstream - 




\ region 


Ji |\|f| S 


1 J ' 











50 



x((im) 



100 



150 



FIG. 1: (Color online) A condensed beam with incident cur- 
rent Ji can populate a stationary scattering state. The solid 
line shows its longitudinal density n(x) (in units of the equi- 
librium density no). In the downstream-region, ip tends to 
a plane wave with transmitted current Jt- The dashed line 
displays the scattering potential V(x) in units of the chemical 
potential u. 



Eq. © exhibits two constant-density solutions: a low- 
density (supersonic) one and a high density (subsonic) 
one, where the transport is respectively dominated by 
the kinetic energy or by mutual interaction of the atoms. 
Both solutions are plane waves of the form ip{x) = Ae lkx , 
but with different wavevectors k and particle densities 
A 2 . As we are considering rather small condensate densi- 
ties and large velocities in the waveguide, the supersonic 
solution will be the relevant one in the context of this 
paper. 

We now assume the presence of a disorder potential 
V(x) in the waveguide which is finite between x = 
and x = L and vanishes elsewhere. In this case, a BEC 
that is injected into the initially condensate-free disorder 
region from the upstream side (i.e. at x < 0) does in 
general not freely propagate to the downstream region 
(at x > L), but undergoes a scattering process. In this 
paper we shall compute transport properties of a system 
where a monochromatic beam of condensate with well- 
defined current Ji is injected into the disorder region (see 
Fig-0. This means, we consider the propagation process 
in terms of a so called fixed input problem |22l |24| . 

Our purpose is now to compute transmission coef- 
ficients for the condensate transport through the dis- 
ordered region. Furthermore, we shall investigate to 
which extent it is possible to populate stationary scat- 
tering states, i. e. stationary solutions ^(cc,t) = 
if)(x) exp(— i/j,t/H) satisfying the outgoing boundary con- 
dition ip(x) = y/noe lkx (with k > 0) for x — » +oo where 
np is the density associated with the supersonic solution 
|29|. This question can be addressed by integrating the 
time-dependent Gross-Pitaevskii equation in presence 
of a source term that is localized in the upstream region 
and emits monochromatic matter waves. Such a source 
models the coupling of the waveguide to a reservoir of 
condensate from which matter waves are injected into 



the guide. It has been demonstrated in that this ap- 
proach is particularly well suited to compute transmis- 
sions for fixed input problems. Additionally, it allows to 
determine for a given potential V(x) whether an incident 
monochromatic beam populates a stationary scattering 
state or not. 

Hence, we consider now the modified Gross-Pitaevskii 
equation with a source which is localized at the position 
xq in the upstream region, 



dip(x, t) 



V(x) 



h 2 d 2 

2m dx 2 
+Sq exp(— i[itjK) S(x 



g\^(x,t)\ 2 



(4) 



So is the source amplitude which determines the emitted 
current. To understand the functionality of the source 
term, it is instructive to consider first solutions of Eq. (QJ 
in absence of the potential V(x). In this case there exists 
plane wave solutions with constant density n. To demon- 
strate this, we switch to the Fourier space, where Eq. Q 
takes the form (for constant n) 



i d 
hdt 



2 JI 



h 2 q 
2m 



gn 4>(q,t) = S e~ 



iqx Q 



-ifit/h 



This equation admits a solution of the form 

So e- iqxo . 



fi — gn — h 2 q 2 /(2m) 



(5) 



(G) 



By transforming back to real space, we find that the 
source emits in both directions the monochromatic wave 
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(7) 



In Eq. J7J k is self consistently defined by (hk) 2 = 
2m [n — g\So\ 2 m 2 /(hrk 2 )] . The current emitted by the 
source can be calculated by evaluating the quantum me- 
chanical current operator. We find Ji — ±|S'o| 2 r7i/(?i 3 fc); 
("+" for x > x , "-" for x < x ). 

We now return to the general case V(x) ^ 0. In order 
to perform the numerical integration, the wavefunction 
ip(x, t) is expanded on a finite lattice and is propagated in 
real time domain. As we are dealing with an open system, 
artificial backscattering at the boundaries of the lattice 
has to be avoided. For that purpose we impose absorbing 
boundary conditions that are well suited for transport 
problems |30| and can be generalized to account for weak 
or moderate nonlinearities [3lj . 

As in real experiments we choose as initial condition 
ip(x,t = 0) = 0. In order to compute the condensate 
wavefunction we numerically integrate ip(x,t) in Eq. 
while adiabatically tuning the source amplitude So from 
up to a given maximum value that corresponds to 
a desired incident current Ji. This approach simulates 
a realistic propagation process, where a coherent Bose- 
Einstein condensate beam with chemical potential fi is 
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injected into the initially empty waveguide from a reser- 
voir. For comparatively weak nonlinearities a stationary 
scattering state of the form ip(x, t) = if)(x) e~ %ixt / fl , which 
corresponds to a supersonic solution in the downstream 
region, is generally obtained from the numerical prop- 
agation. This stationary wavefunction fulfills the time- 
independent Gross-Pitaevskii equation 



(J, tp(x) 



2m dx 2 



■g\Mz)\' 



ip(x) 



(8) 



In contrast to the case of the linear Schrodinger equa- 
tion the transmission coefficient cannot be computed by 
simply decomposing the upstream wavefunction into an 
incident and reflected part because the superposition 
principle is not valid in presence of the nonlinear term. 
Such a decomposition is only possible in the limit of 
small interaction strengths or small back- reflections [32j . 
However, our numerical approach permits nevertheless a 
straightforward access to the transmission coefficient also 
in the nonlinear case. The latter is evaluated by the ra- 
tio of the current Jt in presence of the potential V(x) 
(i.e., the transmitted current) to the current J, obtained 
in absence of V(x) (the incident current that is emitted 
by the source). This approach provides a natural exten- 
sion of the usual definition of transmission coefficients in 
quantum mechanics to the nonlinear case |15| . 

In the nonlinear regime, due to dynamical instabili- 
ties the wavefunction ij)(x, t) does not always converge 
towards a stationary state but can remain time depen- 
dent (cf Sec. 0). In that case, the downstream current 
is no longer constant and therefore the transmission be- 
comes a function of time. In this case we simulate the 
propagation process over a long period r (ideally t —> oo) 
and characterize the transport properties of the guide by 
means of the time-averaged transmission 



T= lim - 

T — >oo T 



t + T 



T(t')dt' 



(*>0) 



(9) 



This choice of working with the mean value T is inspired 
by common experimental setups: the number of con- 
densed atoms Ma reaching the downstream region during 
the time r is Ma — T T J{. This number of atoms can 
be determined experimentally, e.g., by use of absorption 
spectroscopy. 



III. A SIMPLE MODEL OF DISORDER 

In order to compute transport properties through dis- 
ordered regions in magnetic waveguides, it is necessary 
to introduce an appropriate model for the static ran- 
dom magnetic potential along the center of the waveg- 
uide. We first briefly recall the basic principle to gen- 
erate elongated magnetic waveguides for cold atoms or 
condensates. A typical setup that is commonly imple- 
mented on atomic chips is the so called side wire guide 
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FIG. 2: (Color online) Main building block to create a mag- 
netic waveguide on a chip: a current flowing in a microfabri- 
cated wire and a perpendicular bias field form an elongated 
microtrap. Imperfections in the wire force the current to fol- 
low a weakly meandering path and generate therefore a mag- 
netic disorder potential along the center of the guide. 

|33j . As sketched in Fig. 0a circular magnetic field Bo, 
created by an electric current I that flows along a straight 
microfabricated quasi 2-dimensional wire and a homoge- 
neous bias field Bj^ form a minimum of the magnetic 
field parallel to the wire at distance h. An offset field 
B|| applied parallel to the wire reduces losses induced by 
spin flip processes near the magnetic field minimum. For 
a spatially homogeneous current density in an idealized 
wire, the magnetic waveguide is perfectly uniform along 
its longitudinal axis. In reality however, inhomogeneities 
in the current density inside the wire have to be taken 
into account. Such deviations from an homogeneous cur- 
rent flow can be induced by shape fluctuations of the 
wire or impurities inside the metal. These imperfections 
cause a magnetic field roughness along the center of the 
waveguide that acts as an additional potential and pre- 
vents perfect transmission of condensate beam through 
the guide [34| . This additional magnetic field component 
increases as the distance to the chip surface diminishes 
and is expected to reduce the transmission noticeably. 

In the following we consider a steady state current den- 
sity j(r) flowing in a thin quasi 2-dimensional metallic 
wire. Due to the wire imperfections the current density 
varies with the position r. We decompose j(r) into a 
large constant component jo flowing parallel to the wire 
and a small component <5j(r) 



j(r) = jo + <5j(r). 



(10) 



At the center of the waveguide the circular magnetic field 
Bo that is generated by jo cancels with the bias field B^. 
Hence, the total magnetic field along the center of the 
waveguide is given by 



B(x,0,h) 

where <5B = 6B x e x 4 
the Biot-Savart law 

47T 



S|| e x + SB(x,0,h), (11) 
SB z e z is computed from 
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FIG. 3: (Color online) Partitioning of the wire in equal blocks 
of length I, width w, and thickness A z . For each block we 
compute an average current density j". The current compo- 
nent parallel to the ^-direction is at the origin of the magnetic 
disorder potential along the center of the waveguide. 



The effective potential for the atoms is proportional to 
the modulus of the magnetic field 



B| = J(B ]l +SB x 



SB* 



SB 2 ,. 



(13) 



As <5j is supposed to be small, we keep only terms of first 
order in <5j. This yields the simple result 



|B| = B 
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(14) 



Hence, within the approximation of small current fluc- 
tuations, the disorder potential along the center of the 
waveguide is given by 



V(x)=/ j l b SB x (x,0,h). 



(15) 



We now consider a quasi 2-dimensional wire of length 
L in ^-direction and width w in y-direction. A proper de- 
scription of the current density j (r) in the metallic wire 
would require an accurate microscopic model for struc- 
tural dislocations of the wire as well as its impurities 
|34Ll35| . In the present work we adopt a more simple and 
phenomenological approach, which is valid if the length 
scale I on which j(r) typically fluctuates is much smaller 
than the height h of the waveguide. To this end, we di- 
vide the wire into Af equal blocks of length I width w and 
thickness A z (see Fig. [21 and Fig- EJ- F° r each block of 
volume V we compute the average current density 

Az/2 v l w/2 



r = - / dz J dx J dyj(r). (16) 

-Az/2 (f-l)i -w/2 

(The index v = X...M labels the blocks and the corre- 
sponding mean current densities j"). The total electric 
current along the wire is given by 

A z /2 w/2 

1= J dz J dy j(r) • e x = wA z 'f ■ e x . (17) 

-As/2 -w/2 
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FIG. 4: Numerically computed representative examples for 
disorder realizations at different distances h between the cen- 
ter of the guide and the atomic chip surface. The panels 
show the transition from weak to strong disorder, with de- 
creasing distance h. The disorder potential is given in units 
of [io(j,bwAzCt /(2n) 



Hence, in the usual case of a stationary electric current 
/, the x-component of j" is given by the constant value 
jo of Eq. HI U|l for all i/, and we have 



j" = Joe* + 5j»e v + 8f z .e z 



(18) 



The thickness A z of the wire is assumed to be much 
smaller than all other relevant length scales. Wc therefore 
assume \8j%\ <S \Sjy\ and neglect the contribution of Sj z 
to the disorder potential in the following. This yields 



M 



(19) 



where the magnetic field contribution of the v th block at 
the center of the waveguide is computed from the Biot- 
Savart law according to 
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2tt 



arctan 



[(x — x / 

wu/(2h) 
^u 2 + h 2 + w 2 /4 



,/2l3/2 



(20) 
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- x — vl 



Within the discretization procedure just described, 
we can introduce disorder by assuming Sjy as a 
random variable, uniformly distributed in the inter- 
val [-(3cr 2 / 2 1/2 > (3<r 2 / 2 1/2 ]- This assumption cor- 
responds to a zero average transverse current density 
((Sjy) — 0) with a delta correlation 



{Sjy Sjy ) = — $V,V 



(21) 



and allows to reach a well defined regime in the limit 
I — > 0. In that limit Sjy is replaced by an x-dependent 
quantity Sj v (x) verifying jjjfj 



(8j v (x) Sj y (x')) = a 2 S(x - x') 



(22) 
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Here, the parameter a fixes a scale for the typical devia- 
tion of the current density from an homogeneous current 
flow. Since the fluctuations in Sj y are certainly propor- 
tional to jo = I / (wA z ), we can write a = . Here I* 
is a characteristic length depending on the properties of 
the metallic wire, which can in principle be found from 
experimental investigations. 

Due to the convolution procedure in Eq. 1)20(1 . the short 
range disorder in the electric wire induces a smoothly 
varying potential V(x) along the guide. This is clearly 
visible in Fig. 0] which shows the disorder potentials that 
result from three numerically generated sets of current 
densities j" , at three different heights h of the waveguide. 
The disorder potential is smoother for large distances h, 
and becomes more rough (and its typical intensity in- 
creases) as h diminishes. 

One has (V(x)) = and it is appropriate to charac- 
terize the random potential by studying the correlation 
function 
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FIG. 5: (Color online) Numerically computed correlation 
functions for different distances h between the center of the 
guide and the atomic chip surface (solid lines). Dashed lines: 
Fit to a Lorentzian curve. 



C{x-x') = (V(x) V(x')). 



(23) 



In Fig.[S]we show results for C(x-x') at different heights 
h. The correlation function is computed numerically by 
averaging over a large number of different disorder real- 
izations. We find that it can be fitted with good accuracy 
by a Lorentzian curve 



C(x-x') 



lie 



11 + (x - x 1 ) 



IY2 



(24) 



This allows to extract the correlation length l c and to 
establish an empirical relation between the height h and 
l c . In the regime where the width w of the wire and the 
discretization length I are of the same order, we find that 
the correlation length depends linearly on the distance 
between wire and waveguide, l c — ah, with a propor- 
tionality constant s that varies between 1 and 2. For the 
experimentally relevant case of w = 4 fim (a wire of this 
size has been realized by the Tubingen group, see Ref. 
[37|]) we find s = 1.2. Theoretically, this result may be 
understood as follows. In the continuous limit I — > (and 
in the idealized case of an infinitely long wire), Eq. I|2(J|) 
takes, in the regime w S> h, the particularly simple form 
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V(x) 
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, hA z 5j y {x') 
(x - x') 2 + h 2 



(25) 



In this case the disorder potential is exactly Lorentz- 
correlated with l c — 2 h and 7 = 7r~ 1 (poLi B A z o~/2) 2 . 
In the opposite regime w -C h, the correlation function 
C(x—x') cannot be computed analytically, but its Fourier 



transform, C q = J 
lated: One obtains 



+00 



exp(— i q x)C(x)dx, can be calcu- 



^0 fi B wA z a 
2tt 



qKx{qh) 



(26) 



where K\ is the modified Bessel function of the first kind 
[38| . C q as given in Eq. l|26|l is not very different from the 
Fourier transform of a Lorentzian (a decreasing exponen- 
tial), and this is the reason why C(x — x') can be fitted 
reasonably well by a Lorentzian also in the regime /i> w. 
To find a sensible Lorentzian fit, one can for instance try 
to reproduce C(0) and C"(0) obtained from Eq. with 
the parameters 7 and l c of Eq. C(0) = 7/^ and 

C"(0) = -27/^. This leads to 
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+ OC 



2 / t 2 K 2 {t)dt 


+ ft A Ki(t)dt 



1/2 



1.46 . 



(27) 



Thus, again in this limit, we find that l c is of the type l c ~ 
s h. The important outcome of this discussion is that, 
for the continuous model (I — > 0) - in both limits w <C h 
and w 3> h - and also in the numerical realizations of the 
disorder with a finite grid I, we obtain a random potential 
which is Lorentz-like correlated, with a correlation length 
l c that is proportional to the height h of the trap above 
the chip, and with a proportionality constant of the order 
1 — 2. This is confirmed experimentally by the detailed 
studies presented in Ref. [34|. 

The model introduced in this section, where the disor- 
der potential originates from white-noise correlated fluc- 
tuating currents (see Eqs. (|21I22|) ) corresponds physi- 
cally to a "dirty wire model", in the sense that the 
very erratic random current density 11' 21) can be con- 
sidered as originating from the presence of impurities in 
the wire. We note that the white-noise current correla- 
tions lead to a disorder potential whose typical ampli- 
tude (for the continuous limit I — > 0, when h w) varies 
as (V 2 (x)) 1 ^ 2 oc Ih~ 3 / 2 , different from the experimen- 
tal finding I h~ 2 2 of Kraft et dl. [3!|- In contrast, the 
model of "clean wire with corru gat ed boundaries" intro- 
duced in |35| and developed in j 3 ll | yields in the case of 
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a white-noise correlated boundary roughness to a depen- 
dence of the form I h~ 5/>2 , in closer agreement with the 
experimental findings of the Tiibingen group published 
in Ref. |3fj. Note however, that the experimental results 
of the Orsay group point to a boundary roughness 
which is not white-noise correlated, and a typical disor- 
der potential which decreases less rapidly than ft, 5 / 2 , as 
found in the present study. Also the correlation function 
issued from the dirty wire model is in better agreement 
with the experimental one determined in Ref. |34| , which 
differs from the one resulting from a wire with a white 
noise disordered bounda ry ( that has a correlation func- 
tion verifying C 9= o = .35]). It thus appears that the 
simple dirty wire model introduced in the present section 
allows to construct a disordered potential V(x) that cap- 
tures most of the characteristics of the micro fabricated 
magnetic guides. 



rewrite Eq. I|30|) in a dimensionless form: 



i dV /2 p-1 



p 



1/2 5a; 2 



-- 1 

P 2 



v(x) 



= 0. (31) 



In this expression we made use of the dispersion rela- 
tion Eq. © and expressed J — n^hk/m in terms of the 
downstream density n and of the outgoing wave vector 
k. The quantity H 2 k 2 /(2m) = l-i — griQ is the kinetic en- 
ergy of the outgoing plane wave with equilibrium density 
riQ. £ = h/^/2mnog is the condensate healing length. 

In order to find perturbative solutions of Eq. I|31|l for 
v(x) -C 1, we insert the ansatz p{x) = 1 + 5p[x) into 
Eq. and keep only terms that are linear in Sp(x): 

d 2 

-7Ti5p{x) + 4k 2 6 P (x) = 2k 2 v(x) , (32) 
ox A 

where 



IV. WEAK DISORDER 

In this section we investigate the regime of weak 
disorder potentials and derive simple relations between 
the condensate density and the disorder potential V(x). 
Weak disorder means in this context that the propaga- 
tion of the condensate is only marginally affected by the 
scattering region. This implies that the kinetic energy 
per particle must be much larger than the typical inten- 
sity of the disorder potential (which can be estimated 
for instance by the standard deviation (V 2 (x)) 1 ^ 2 ). We 
shall argue below that a secondary criterion is necessary 
to characterize this regime, namely that the length of 
the disordered region is small compared to the charac- 
teristic length scale Lj, (to be defined below) typical for 
the decrease of the transmission. 

First, we rewrite the Gross-Pitaevskii equation in 
the well known form of hydrodynamic equations 



d 

dl 1 



d_ 

dx 



(nu) 



(28) 



and 



du 



d_ 

dx 



d 2 n x l 2 mu 2 



2mn 1 / 2 dx 2 



V(x) — gn 



(29) 

where u is the condensate velocity. In the case of a sta- 
tionary state we have <9 t n = and dtu = from which it 
follows that the current J = nu is constant. Integration 
of Eq. (£211 then yields 



,i = V(x) 



gn 



J 2 



h 2 d 2 ^' 2 



2mn 2 2mn 1 / 2 dx 2 



(30) 



This is the time-independent Gross-Pitaevskii equation 
for a current carrying scattering state. In the down- 
stream region an outgoing plane wave ip(x) = y/noe lkx is 
expected. The equilibrium density no coincides with the 
supersonic solution of the dispersion relation Q . Defin- 
ing p(x) = n(x)/no and v(x) = 2mV(x) / {h 2 k 2 ) one may 



k = k\ / 1 — 



2£ 2 fc 2 ' 



(33) 



The solution of Eq. I|32|) in presence of the downstream 
boundary conditions 5p(L) = 0, Sp'(L) = (flat down- 
stream density) is |4l| 



5p(x) 



K 



2k(x' — x) v(x') dx' , 



8p'{x) = -2k 2 / cos 2k{x'-x) v{x')dx' . (34) 



This implies that the density profile in the upstream re- 
gion (x < 0) deduced from the linearized Eq. l|3^|) is of 
the form n(x) = uq [1 + Sp(x)] with 



Sp(x) = Sp cos(2kx + 9) 



(35) 



The amplitude Sp and the phase factor 8 in Eq. I|35(l are 
determined by the disorder potential V(x) via Eq. J22J. 
The modified wave number k fixes the period of the den- 
sity oscillations. 

As we are obviously in the regime of small back- 
reflections we adopt the method of Ref. [12 to determine 
the transmission coefficient in an approximative way. To 
this end we make the ansatz n(x) = |^ inc (a;) + ip„ { {x)\ 2 
with 



= a exp{inx} , 
ip IB i{x) = b exp{i(ftx + 9)} 



(36) 



Comparing the corresponding density profile with ||5SJ|, 
one obtains the following expressions for the amplitudes 
a and b [40| : 



1 



"o 
b 2 
no 



Sp 2 



-ST 



0(Sp 4 ) 



(37) 
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It was pointed out in Ref. (33 , and numerically confirmed 
for single and double barrier potentials that ip Tcl 

can be approximately identified with the reflected com- 
ponent of the condensate in the case of almost perfect 
transmission. This corresponds to a reflexion coefficient 
R = b 2 /a 2 = j S~p 2 + 0(5~p 4 ) and to a transmission coef- 
ficient which can be expressed (using Eq. I|35|)l as 



T = 1- 



;Sp 2 



4 K 



[Sp'(0)} 2 . (38) 



In this final expression Sp(Q) and 5p'(0) are related to 
the disordered potential by means of Eq. Therefore, 
determining the transmission T for a given potential V(x) 
amounts to compute the integrals Eq. pijl. 

As shown in AppendixEJthe above procedure allows to 
determine the disorder average (T) from the knowledge 
of the correlation function of the disorder potential. For 
the relevant case of a Lorentzian correlation (of the form 
pi)ll we obtain 



(T) = 1 - 



L 



where 



L d = 



h^ 2 

iTjm 2 



(39) 



(40) 



is the characteristic length scale for the decay of the 
transmission. We recall here that the above analysis is 
only valid in the regime <5p <C 1, i.e., the linear decrease 
of (T) in Eq. |J55J| is only valid for L <C L d . Thus, we 
have to refine our definition of weak disorder: not only 
the intensity of the potential should be small, but also 
the length of the disordered region should not exceed the 
value Ld- 

As we see from expression l|40|) . the effect of the atom- 
atom interaction is entirely contained within the modified 
wave number k (Eq. 1)33(1 ^ which describes the period of 
the upstream density oscillations. For repulsive atom- 
atom interactions, we have k < fc, which implies that 
the mean transmission is reduced compared to the non- 
interacting case. This behavior is indeed well confirmed 
by numerical computations based on the approach pre- 
sented in Sec. [3 This interaction-induced decrease of 
the transmission was already observed in Ref. ^ij^ and 
interpreted as a lack of kinetic energy compared to the 
interaction-free case. 

In the limit of very small correlation lengths, i.e. 
kI c <C 1, the disorder potential can be approximated 
by a white noise potential with correlation function 
C(x — x') ~ "fnS(x — x'). Considering the noninteracting 
case (k = k) we recover in this regime the well-known ex- 
pression Ld — Li oc = (h 4 k 2 )/ (irm 2 ^) for the localization 
length of ^-correlated disorder potentials (see, e.g., Ref. 

Ef). 

The opposite limit kI c ^> 1 can be considered as the 
semiclassical regime, where the de Broglie wavelength 
A = 2ir/k of the condensate is much smaller than the 



correlation length l c of the disorder potential. In this 
regime, the length scale Ld is dominated by the exponen- 
tial prefactor exp(2«;/ c ), and the deviations from perfect 
transmission (T) = 1 vanish exponentially fast with in- 
creasing ratio kI c . The semiclassical condition kI c 3> 1 
furthermore allows to derive a simple analytical expres- 
sion for the density n(x) throughout the scattering re- 
gion. We start from the zeroth-order solution n(x) = uq 
valid for V = 0. Then, for given /i and J the density no 
can be obtained by iteratively solving the selfconsistent 
equation (strictly equivalent to J3J)) 



n 



J [p - gn ] 



(41) 



starting, e.g., with n — J^m/ (2/i). This procedure 
guarantees convergence toward the supersonic solution 
of Eq. ®. 

The natural generalization of Eq. 141fl to the case of 
a small but non- vanishing potential V(x) is obtained by 
using Eq. instead of Eq. This yields 



n = uq 



V{x) V-gn+z^g^- 



p - gn Q 



p - gn 



where the current J was substituted by means of the dis- 
persion relation @. We shall now find approximate so- 
lutions of this selfconsistent equation in the case of weak 
disorder, i.e. \v(x)\ -C 1, where the typical value of V is 
much smaller than /i — gn which is the kinetic energy 
per particle. We emphasize that this does not imply that 
the nonlinear term griQ should be small. 

The zeroth order solution of Eq. I|42|) is simply the 
constant equilibrium density tiq. Resubstituting this con- 
stant solution into the recursive equations yields the first- 
order solution for the condensate density 



(x) 



"0 



y/l-v(x) 



(43) 



Corrections to this first-order expression particularly 
arise from the quantum pressure term 9 ^ ^ ^2\/w- ^ 
can be shown, however, that the latter is suppressed by 
a factor ~ l/(kl c ) 2 as compared to the kinetic energy 
H 2 k 2 /{2m) when n^ x \x) is resubstituted in Eq. (|4*2*)) . In 
the semiclassical regime kl c 3> 1, the quantum pressure 
term becomes negligible, and the expression (|43(l repre- 
sents a very good approximation to the actual density of 
the condensate in the scattering region. We show in Ap- 
pendix[B]that the result (|43|l can be derived in a way that 
is directly analogous to the semiclassical WKB approach. 

This result is illustrated in Fig. whose lower panel 
shows a random potential generated with the method 
presented in the previous section. In the upper panel we 
compare the result of Eq. Q43|) with an exact, i.e. numer- 
ically computed solution of the Gross-Pitaevskii equa- 
tion. Excellent agreement between the first order solu- 
tion and the exact solution is found. We note here that 
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x(p.m) 

FIG. 6: (Color online) The upper panel displays a comparison 
of the first order solution (WKB, Eq. I43H ) with a numerically 
computed solution (QM) of the Gross-Pitaevskii equation for 
a weak disorder potential V(x) (shown in the lower panel). 
The correlation length is l c — 30 /im, the wavelength is A = 
3 fim. The ratio between interaction and kinetic energy in the 
incident beam is Ei nt /Ekin = 1/10. 



it is quite natural to find that the density profile mir- 
rors the potential because we are dealing with current 
carrying states: the relation (|30f) (in absence of quantum 
correction) /i = rau 2 /2 + V(x) + gn predicts that the 
condensate velocity becomes minimal close to the max- 
ima of the disorder potential. It follows then from the 
continuity equation J = nu = const that the density n 
assumes its maxima when the velocity becomes minimal. 

It is instructive to realize that classically forbidden 
back-reflections can be taken into account by inserting 
the ansatz n(x) — n^(x) + Sn(x) into Eq. (|30|l and lin- 
earizing the resulting equation for small 5n{x)/riQ. To 
the lowest non vanishing order in v, we again obtain the 
result 139fl for the mean transmission. 

Finally we consider experimental realizations of waveg- 
uides on atom chips. Typical distances h between the 
chip surface and the guide are in the range 20 - 100 /im. 
Typical disorder correlation lengths are of the same order 
as h. In recent transport experiments 01 the velocity of 
propagating 87 Rb condensates is of the order of a few 
millimeters per second, resulting to a mean wavelength 
of a few micrometers. This corresponds to the regime 
kI c ^> 1 with - from (|4"U)l - a very large value of L d . 
Hence the regime of weak disorder is presently the most 
relevant one: the kinetic energy is much larger than the 
typical intensity of the disordered potential and L e i is 
large compared to the typical length of the disordered 
region: one thus expects almost perfect transmission. 



V. MODERATE AND STRONG DISORDER 

In Sec. IIVI we focused on weak disorder potentials, in 
the limit of small reflection. The analysis was done in the 
regime fx ^> (V 2 (x)) 1 ^ 2 and L <§; Ld- In the present sec- 



tion we still partially fulfill the first of these inequalities, 
but drop the second one. We will see that the behavior 
of the system is quite different, ranging from a regime of 
localization (in the limit of weak interaction), to a time 
dependent behavior for larger interaction, with a power 
law decay of the time-averaged transmission. 

First, we shall discuss some elementary differences be- 
tween the scattering problem in linear quantum mechan- 
ics and the nonlinear Gross-Pitaevskii equation. In linear 
quantum mechanics (g = 0) one finds for any scattering 
potential a unique stationary scattering state that is dy- 
namically stable, and the associated transmission coeffi- 
cient T relates the constant incident current Ji one-to- 
one with the transmitted current J t . For the nonlinear 
Gross-Pitaevskii equation the transmission T depends on 
the density of the propagated condensate and thereby on 
the current. Additionally, the phenomenon of multista- 
bility may arise. This means that for a given incident 
current Ji two or more scattering states with different 
transmissions can coexist. 

In principle all stationary scattering states that are as- 
sociated with a given incident current Ji can be found by 
integrating the time-independent Gross-Pitaevskii equa- 
tion JSJ) from the downstream to the upstream region. A 
systematic variation of the downstream current Jt allows 
to select the desired states. This procedure, however, 
does not reveal any information about their dynamical 
stability properties, which are crucial for answering the 
question whether an incident condensate beam populates 
a stationary scattering state or not. For instance, in the 
case of coherent condensate transport through a double 
barrier potential, three possible scattering states are ex- 
pected close to the resonances, but only one of them is dy- 
namically stable Here the advantage of integrating 
the time-dependent Gross-Pitaevskii equation becomes 
apparent: If this integration converges to a stationary 
scattering state we know automatically that this state is 
dynamically stable (otherwise small numerical deviations 
would exponentially increase with propagation time) . 

We consider an ensemble of ./V disorder realizations 
with randomly varying sample lengths L that are uni- 
formly distributed between and a maximal sample 
length. For each realization (labeled with index a) we 
numerically compute the time evolution of the wavefunc- 
tion and extract either the time-independent transmis- 
sion T a (if if}(x,t) converges to a stationary state) or the 
time averaged transmission T a (if ip(x,t) remains time- 
dependent). For the sake of definiteness and due to its ex- 
perimental relevance we consider the propagation of con- 
densed 87 Rb atoms (whose scattering length is a sc = 5.77 
nm). Our numerical computations were performed for a 
guide with radial trapping frequency u>± = 2ir x 100 s _1 
(oscillator length a± = 1 /im) . The disorder is generated 
as in the previous section. The regime of strong disorder 
is reached by choosing a rather short distance h = 5 /im 
between the center of the guide and the chip surface, 
which corresponds to a correlation length l c = 6 fim. In 
order to avoid excitations of the condensate into higher 
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transversal modes we adjust the standard deviation of 
the potential (which is a measure of the mean potential 
height) to {V 2 {x)) 1 / 2 ~ 0.12 hujj_. In allthe following nu- 
merical calculations we consider an incident monochro- 
matic beam with current Ji — 10 3 atoms per second and 
wavelength A = 10 fj.rn. Then the chemical potential is 
ii = 0.25 Hu>± (in the linear case the chemical potential 
takes the slightly different value /i = 0.23 hui±). 

It is instructive to focus first on the linear case (g — 
0) that has already been extensively investigated in the 
context of localization theory [22L . In the localized 
regime the transmission decays exponentially with the 
system length L, i.e., (T) = exp(— L j Li oc ) where L; oc is 
the so called localization length. The points in the upper 
panel of Fig. mark for each disorder realization the 
associated transmission T a (L) as a function of the sample 
length L. To extract from these data a characteristic 
scaling law for the L-dependence of the transmission we 
divide L into equal intervals of length AL <§; L. We then 
compute the mean transmission at sample length L by 
summing up all the values T a corresponding to a sample 
length lying in the interval of width AL centered at L: 



(T)a(L) 



Nl is the number of samples in the interval under con- 
sideration. 

The step function in Fig. [7| shows the decrease of (T) a 
for 30000 disorder realizations and AL = 50 (for the 
sake of clarity we show only 2000 points in the plot). 
In the context of localization theory it is convenient to 
investigate scaling laws by means of the geometrically 
averaged transmission 



,<MT)) 



(ln(T)) 



Lj2ln(T a (L')). 



(45) 
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FIG. 7: (Color online) Transmission through a disor- 
dered sample as a function of sample length L for the non- 
interacting case. Each point corresponds to a different real- 
ization of the disordered potential. Upper panel: arithmeti- 
cally averaged transmission (blue staircase function). Lower 
panel: The geometric averaged transmission (blue staircase 
function) decreases exponentially with L, as revealed by the 
fit with Lioc — 586 /im (straight black line). The arrows mark 
the logarithmic standard deviation. 



because, contrary to (T) a , the average (ln(T)) is a self- 
averaging quantity of the system |42t |43j. The lower 
panel of Fig. [7| shows (T) g which follows clearly an expo- 
nential law. This is a clear evidence for the appearance 
of localization. We can extract the localization length 
that is here L/ oc = 586 /im. We note the wide spread 
of the data points around their average. This spread is 
quantified by the logarithmic standard deviation 



Aln(T) 



Nr 



£[ln(T a )-(ln(T)) 



(46) 



which is shown as arrows in the lower panel of in Fig. \7\ 
We find an almost linear increase of Aln(T) with the 
sample length. 

Is the conventional localization scenario, with the 
characteristic exponential decrease of the transmission 
[l6l , still valid in the case of interacting atoms? To 
address this question we now calculate the transport in 
presence of a moderate nonlinearity where the ratio of 



interaction and kinetic energy in the incident beam is 
Eint/Ekin — 1/10. Contrarily to the linear case, time- 
dependent behavior becomes now a dominant feature as 
shown in Fig. El We find that dynamical stable scattering 
states (black dots in Fig.[8J) are only populated for sample 
lengths that are smaller than a critical length L* which is 
here of order of 125 /xm. For samples with length L >, L* 
we find a crossover region where time-dependent behav- 
ior sets in and convergency to a stationary state is only 
achieved for a certain fraction of disorder samples. ip{x, t) 
remains time-dependent (orange crosses) for all samples 
when we reach the regime where L is notedly larger than 
L* . In the time-dependent case the data points display 
the time averaged transmissions T a 

In order to extract a scaling law from our data, we com- 
pute the ensemble-averaged transmission (in the time- 
dependent cases T a in Eqs. 144I46|I is replaced by T a ). 
We find that the geometric averaged transmission (T) g 
(step function in Fig.[8J) decreases inversely with the sam- 
ple length L and is well approximated by the algebraic 




L(nm) 

FIG. 8: (Color online) Transmissions through an ensemble of 
disorder realizations for a moderate nonlinearity Ei nt / Ekin = 
1/10; The characteristics of the incident beam are given in the 
main text. A transition from a time-independent to a time- 
dependent regime is observed at critical length L* ~ 125 \i 
m. The black dots represent the transmissions in the time- 
independent regime and the (orange online) crosses are the 
time averaged transmissions in the time dependent regime. 
The staircase function is the geometrically averaged trans- 
mission. It is well approximated by the algebraic scaling law 
Lo/(L + Lo) (smooth solid line) with Lq = 287/im. 

function (smooth line in Fig. [SJ 



with the decay length Lq. Such a scaling law is character- 
istic for transport in systems with loss of phase coherence 
between the single scattering events. Indeed, if one con- 
siders a series of successive scatterers and calculates the 
transmission by neglecting all interference effects one de- 
rives exactly the scaling law of Eq. l(4"T|l [ill Iit3 | . Such 
an ohmic behaviour is observed for electron transport 
through mcsoscopic metal structures in the limit of small 
dephasing lengths (4^, |4g . 

Another striking feature is the distribution of the data 
points in Fig. |SJ Contrarily to the linear case, this dis- 
tribution is now clearly restricted to the neighborhood 
of the average transmission and the standard deviation 
AT(L) decreases for long sample lengths L. Hence, in the 
regime of large lengths one expects to find the T a 's in a 
narrow interval centered around the averaged transmis- 
sion. Loosely speaking, T a becomes more or less sample 
independent. For the sake of completeness we mention 
that ideally the time averages T a should be computed for 
an infinitely long period. Of course this cannot be done 
numerically, but we verified that the averaged transmis- 
sion and the standard deviation do not change if we in- 
crease in (01 the averaging time window from r to 2r and 
3t. 

The above presented computations demonstrate that 
even a moderate nonlinearity leads to a dramatic change 
of the transmission properties. In particular, the usual 
interpretation of the transmission behavior in terms of 




FIG. 9: (Color online) Transition from time-independent to 
time-dependent behavior in presence of a very weak nonlinear- 
ity (Eint I Ekin = 1/100). The staircase function in the upper 
panel shows the geometric average of both time-independent 
and time-dependent transmissions (black dots and orange 
crosses, repsectively) . The lower panels display separate av- 
erages over the time-independent transmissions (lower left 
panel) and the time-dependent transmissions (lower right 
panel) and show the best exponential (dashed line, colored 
red online) as well as algebraic (solid line, colored red on- 
line) fits to the data. Clearly, the time-independent transmis- 
sions decrease exponentially with L (with localization length 
Lioc = 439 ^tm) while the averaged transmissions in the time- 
dependent case decay according to Ohm's law. 



localization is put at question in the case of interact- 
ing particles. In order to obtain deeper insight into that 
matter, we redo the above computation with a very weak 
nonlinearity, such that E int /Eki n = 1/100. Fig. |5] shows 
that for this case the crossover from time-independent 
to time-dependent behavior is shifted to larger sample 
lengths (L > 600 fim) . This indicates the existence of a 
critical nonlinearity beyond which the system exhibits 
time dependence. Indeed, preliminary studies show that 
for each disorder sample length L there is a critical value 
g* above which no stationary scattering state can be pop- 
ulated, or, equivalently, for each strength of interaction 
<7, there is a critical disorder length L* above which the 
flow is time dependent. We find that L* decreases with 
increasing nonlinearity. This is the reason why station- 
ary states can be populated in short but not in long dis- 
order regions and why the crossover to time-dependent 
dynamics in Fig. |5| is displaced to larger sample lengths 
L compared with Fig. |S] 

From Fig. [5] we also see that the Ohmic decrease of 
the transmission is intimately connected to the occur- 
rence of time-dependent dynamics. This is clearly indi- 
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cated by the fact that the time-independent data points 
can be accurately fitted by an exponential law, and not 
by an algebraic one [^. We infer from this observa- 
tion that as long as stationary states are populated the 
system follows the conventional scenario of localization 
even in the presence of repulsive atom-atom interactions, 
with a smaller localization length than in the interaction- 
free case. This scenario seems to break down as soon 
as the scattering process of the condensate becomes in- 
trinsically time-dependent. We tentatively attribute this 
phenomenon to the fact that the definition of the mean 
transmission involves a time average over the propaga- 
tion process (see Eq. (JSJ). Therefore, information about 
the phase coherence, which is in principle preserved by 
the time-dependent Gross-Pitaevskii equation, becomes 
lost in the time-averaging procedure. 



VI. CONCLUSIONS 

In this article we have presented a study of transport of 
Bose-Einstein condensates in presence of disorder. We in- 
troduced a one dimensional model for the disorder poten- 
tial in the case of a condensate that propagates through 
a magnetic waveguide over a microchip. We assumed for 
this model that the transverse current density in the mi- 
crofabricated wire exhibits a white-noise correlation. We 
could show that this yields a disorder potential that is 
Lorentz-like correlated along the axis of the waveguide. 

In the regime of weak disorder, a perturbative ap- 
proach allowed to estimate the deviations from perfect 
transmission. We found that on length scales much 
smaller than a characteristic length L c i - which is de- 
termined by the correlation length of the disorder po- 
tential and the healing length of the condensate - the 
transmission decreases linearly with the length L of the 
disorder region. The presence of a repulsive atom-atom 
interaction diminishes the transmission compared to the 
interaction-free case. Furthermore, in the limit of large 
correlation lengths (l c 3> A) we could identify a semi- 
classical regime where the backscattering is exponentially 
suppressed and where the condensate density mirrors the 
shape of the disorder potential V(x). 

The numerical approach presented in Sec. ITT1 provides 
an access also to the regime of moderate and strong disor- 
der potentials and allows to simulate a realistic transport 
process. In the case of noninteracting atoms we find clear 
evidences of the appearance of localization. In presence 
of interaction, time-dependence of the transmissions be- 
comes a dominant feature of the system. We find that 
stationary scattering states can then only be populated 
in waveguides with rather short disorder regions, whereas 
the condensate exhibits a strongly time-dependent dy- 
namics if we consider large sample lengths. Our nu- 
merical calculations show that the critical length L* , at 
which the crossover between the two different regimes oc- 
curs, is shifted towards shorter sample lengths when the 
strength of the interaction is increased. It remains an 



open problem to determine this critical length L* analyt- 
ically from the system-specific parameters which are the 
incident current, the kinetic energy, the average height 
of the disorder potential, and the associated correlation 
length. 

Our numerical study was restricted to Lorentz- 
correlated disorder in the atom-chip context. We ex- 
pect, however, no significant differences for other types of 
disorder potentials, such as speckle fields or point scat- 
terers. Indeed, preliminary studies on the transport of 
BEC in presence of randomly placed <5-like barriers re- 
veal qualitatively the same phenomenology: A regime of 
time-dependent scattering sets in beyond a critical inter- 
action strength (or sample length), and the transmission 
decreases according to an Ohm-like law rather than to 
an exponential one. An important aspect that remains 
unexplored, on the other hand, is the depletion of the 
condensate, due to atom-atom scattering events, and the 
appearance of a thermal cloud that propagates together 
with the condensate. This issue should be rather relevant 
in the regime of time-dependent scattering, but cannot 
be studied with our present approach which is based on 
the mean-field description of the condensate. Since the 
interaction with such a thermal cloud will lead to an addi- 
tional cause for incoherent transport, we expect no quali- 
tative change as far as the Ohmic power-law decay of the 
transmission is concerned. 

The present work opens new perspectives for the study 
of transport in phase coherent systems. In particular the 
repulsion between the atoms leads to a behavior different 
from the one expected in the non interacting and attrac- 
tive cases. Throughout this work we have considered 
realistic values of the parameters (intensity and correla- 
tion length of the potential, distance from the guide to 
the microchip, incident current of the beam) describing 
a BEC of 87 Rb in a waveguide and hope to motivate ex- 
perimental studies testing the results presented in this 
work. 
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APPENDIX A 



In this appendix we derive a relation between the mean 
transmission (T) and the correlation function of the dis- 
order potential in the weak disorder limit. Taking the 
mean value of Eq. (3HJ gives 



(T) = l-\l([5 P (0)] 2 ) + ^- I ([Sp'(0)] 2 )} . (All 



Therefore the problem of computing (T) reduces to the 
calculation of the averaged values 



<[M0)] 2 } = ^E^(°)] 2 ' 

i=l 
1 N 

([Sp(0)] 2 ) = -J2M(0)} 2 - (A2) 



Eq. H34|) allows us to write these averages as 



an exact analytic evaluation of the integrals (|A3|) is not 
possible. Nevertheless, in the regime L l c and L ^> 
1/k the inner integral can be approximated with high 
accuracy by the real or respectively the imaginary part 
of 



1 L e 



2inx 



dx' 



7T I 2 + (x ~ X') 2 



; 2 «(«-« [e(a?) - Q( x - L)} 



The same integration by parts procedure as in the white 
noise potential yields a result that is only modified by the 
occurrence of the exponential factor exp(— 2nl c ). There- 
fore, in presence of a Lorentzian correlation the mean 
transmission reads 



(T) = 1 - -^e-«. 



APPENDIX B 



(A9) 



([*p(o)r> = - x 

„ . 



([V(0)] 2 ) = 4fc 2 x 

L L 

x J dxcos{2nx) J dx' cos(2kx')(v(x)v(x')) ,(A3) 
o o 

where (v(x)v(x 1 )} is the correlation function of the po- 
tential. We evaluate (T) for the particularly interesting 
cases of a delta-correlated white-noise potential and a 
Lorentzian-correlated disorder potential. In the first case 
the correlation function reads 

( V ( x ) v ( x >)) =1 ^—j nS(x-x'), (A4) 

and the integrals in Eq. (|A3|I can be easily evaluated by 
means of integration by parts. This yields 



2 2?7i 2 77r 



L + c 



(A5) 



where c is a dimensionless constant which is of the order 
(477T771 2 ) / (h 4 k 3 ) . For sample lengths L ^> we can 
neglect c and we keep only terms that scale linearly with 
L. By using this approximation we find with Eq. (|A1|) 
for a white-noise random potential 



(T) = 1 - 



irm 2 jL 
h A n 2 



In the case of a Lorentzian correlation 

2 



{v(x)v(x')) = 7 



/ 2m 



\h 2 k 2 J l 2 c + {x-x<) 21 



(A6) 



(A7) 



In this appendix we demonstrate that the first order 
solution Eq. (|43ll can also be derived from a standard 
WKB - method |49|. Inserting the WKB ansatz tp(x) = 
exp(i/(x)) (where f{x) is a complex-valued function) in 
the time-independent Gross-Pitaevskii Eq. (JS) gives 



~ [if -iff] 



+ V(x) + ge 



Kf-n 



/i. 



(Bl) 



The potential V{x) is supposed to be small compared to 
p — gn , and its correlation length is large compared to 
the de Broglie wavelength. Eq. (|B1|) can be rewritten in 
terms of a selfconsistent equation 

/' = yJ^{»-V{x)-gjU-n)+if». (B2) 

which allows to compute recursively the unknown func- 
tion f{x). Treating V(x) as a small perturbation moti- 
vates to use a plane wave with wavenumber k and density 
no as zeroth-order approximation 



/(°)=fcx-iln(Vn7) 
Hence, by use of Eq 1B2|) . we find 



(B3) 



/«' = ^y/2m(n-V(x)-gn ) = ±P(x), (B4) 

and the first-order approximation of the wavefunction 
reads 



/ X \ 

^ {1) (x) = y/h~o~ exp j i J ^P(x')dx' 



(B5) 



The lower bound a of the integral lies in the downstream 
region, where we assume the flat condensate density no. 
Inserting the second derivative 



/ 



(i)' 



~KT\ V '& = \P'(x) (B6) 
h P(x) n 
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into the selfconsistent Eq. i|B2(l yields the second-order 
approximation 



Hence, the density found with this WKB method is 



CO' = ^L^ l + lh p l{x) / { p [x) y 



h 2 

P(x) P'(x) 



2P(x) 



(B7) 



By integrating we find the second-order approximation 
of the wavefunction 



^ 2 \x) = exp ( i J \p{x')dx' ) . (B8) 



n(x) 



no 



Ha) 

P(x) ^l-v(x) 



(B9) 



where v(x) = V(x)/(fi — gno) = V (x)2m / (ft 2 k 2 ) . Indeed, 
Eq. 1|B9|) coincides exactly with the result (|43J) found in 
Sec.Hvl 
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